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ABSTRACT 

This  paper  presents  the  performance  analysis  of  two 
digital  beam  forming  techniques  used  in  conjunction  with 
a  software  GPS  receiver  to  mitigate  interference  to  GPS 
signals  in  jamming  environment.  The  first  method  is  the 
constrained  minimum  power  (MOP)  method.  The  second 
method  is  the  so-called  self-coherence  restoral  (SCORE) 
method.  Both  experimental  and  simulation  data  are  used 
in  the  study. 

The  study  was  performed  using  experiment  data 
collected  in  an  anechoic  chamber  to  obtain  GPS  and 
interference  signals.  A  two  by  two  GPS  antenna  array 
and  a  four  channel  radio  frequency  front  end  are  used  to 
collect  simulated  GPS  data  generated  using  hardware- 
based  simulator  in  controlled  jamming  environment. 
Three  types  of  interference  signals  are  deployed  in  the 
experiments:  EM  chirp,  binary  phase  shift  key,  and 
broadband.  The  interference  power  levels  used  were  +20, 
+30,  and  +40  dB  above  GPS  signal  power.  A  software 
GPS  receiver  was  used  to  perform  acquisition  of  GPS 
signals  to  evaluate  the  performance  of  the  beam  forming 
algorithms.  The  preliminary  result  showed  that  MOP 
method  can  effectively  mitigate  all  three  types  of 
interference  at  all  power  levels  if  a  single  jammer  is 
present.  Experiments  using  multiple  broadband  jammers 
were  also  analyzed  and  our  results  shown  that  the 
effectiveness  of  the  MOP  diminishes  as  the  jammer  power 
increases  and  ceases  to  function  at  the  +40  dB  level.  The 
SCORE  method  does  not  exhibit  consistent  performance 
for  the  experimental  data.  This  is  consistent  with  our 
simulation  results  which  show  that  for  the  SCORE 
algorithm  to  generate  satisfactory  results,  sufficient 
number  of  antenna  elements  is  necessary  even  if  there  is 
no  interference  source  present.  The  number  of  antenna 
element  is  determined  by  the  number  of  satellites 
available,  as  well  as  the  number  of  interference  source. 
Detailed  analysis  of  the  simulation  and  experimental 
results  are  discussed  in  this  paper. 
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1.  INTRODUCTION 

A  number  of  methods  and  approached  have  been 
proposed  and  reported  to  mitigate  interference  for  the 
GPS  receivers  in  jamming  environment.  Most  of  these 
reports,  however,  are  based  on  simulation  results  (Amin  et 
al,  2003;  Amin  et  al,  to  appear;  Blazquez  et  al,  1999; 
Fante  and  Vacarro,  1998,  2000;  Hatke,  1998;  Moelker  et 
al,  1996;  Zhang  et  al,  2001;  Zoltowski  and  Gecan,  1995). 
Recently,  efforts  are  being  made  to  investigate  the 
experimental  performance  of  the  various  beaming 
forming  techniques  (Liou  et  al.,  2001;  McDonald  et  al, 
2004).  In  this  study,  two  beam  forming  techniques,  the 
constrained  minimum  power  (MOP)  method  and  self¬ 
coherence  restoral  (SCORE)  method,  were  implemented 
to  evaluate  their  claimed  performance  using  experimental 
data  collected  in  controlled  jamming  environment. 
Simulation  studies  were  also  performed  to  validate  the 
algorithms.  A  software  GPS  receiver  is  used  to  perform 
GPS  signal  acquisition  of  the  beam  former  output. 

The  MOP  method  utilizes  the  fact  that  GPS  signals 
are  far  beneath  the  thermal  noise  level.  Minimization  of 
total  GPS  receiver  input  power  while  maintaining  gains 
along  the  directions  of  GPS  satellites  will  suppress  the 
interference  contribution.  The  SCORE  method  is  based 
on  known  repetitive  nature  of  the  GPS  signal  CA  code. 
GPS  signal  samples  separated  by  integer  multiples  of  the 
CA  code  length  while  within  the  same  navigation  data  bit 
are  self-coherent  because  of  the  code  repetition  property. 
Both  methods  have  claimed  pros  and  cons  in  the 
literature.  MOP  method  may  work  effectively  on  a 
variety  of  jamming  sources.  But  it  requires  prior 
knowledge  of  the  satellite  orientation.  The  SCORE 
method  does  not  need  any  satellite  position  information. 
It  is,  however,  ineffective  against  jamming  sources  that 
have  self-coherent  prosperities.  Moreover,  our  study 
showed  that  for  the  SCORE  method,  the  number  of 
antenna  element  is  determined  not  only  by  the  number  of 
interference  sources,  but  also  by  the  number  of  satellites 
available. 

Section  2  of  the  report  will  describe  the  experimental 
setup  that  generated  the  data  for  the  study.  Section  3 
summaries  the  MOP  and  SCORE  methods  implemented 
in  the  study.  The  simulation  results  and  experimental  data 
analysis  will  be  presented  in  Section  3  and  4  respectively. 
Section  5  concludes  the  summer  project  and  discusses 
future  works. 

2.  EXPERIMENTAL  SETUPS 


1  through  4  respectively.  The  distance  d  from  an  antenna 
element  to  the  center  of  the  dome  is  7  cm.  Phase 
calibration  of  the  antenna  elements  has  been  performed 
and  the  results  were  used  in  the  beam  forming  algorithm 
(Liou  et  al,  2002).  The  distance  between  the  elements, 
the  size  of  the  antenna,  and  the  size  of  the  geodesic  dome 
are  not  drawn  to  proportion. 


Fig.  1.  Antenna  elements  layout  and  experimental  setup. 

Eight  GPS  antennas  transmitted  signals  from  a  GPS 
simulator  were  mounted  on  the  dome  to  mimic  a  GPS 
constellation.  Three  extra  antennas  were  placed  nearby  to 
transmit  interference  signals.  The  interference  signal  has 
three  modulation  schemes:  a  bi-phase  shift  key  (BPSK) 
source  with  10k  and  IM  modulation  rates,  a  frequency- 
modulated  (FM)  chirp  source  with  Ik  and  10k  modulation 
rates,  and  a  35  MHz-bandwidth  broadband  random  noise. 
The  power  levels  were  +20,  +30  and  +40  dB  above  the 
GPS  signals  in  the  El  band.  A  sky  plot  of  the  GPS 
transmitters  and  the  interference  sources  placements  are 
shown  in  Figure  2.  The  interference  sources  are 
represented  by  Jl,  J2,  and  J3,  while  the  numbers  4,  5,  6, 
10,  13,  24,  26,  and  30  represent  the  GPS  transmitters, 
with  the  numbers  representing  GPS  satellite  ID.  The 
azimuth  angle  ^  is  measured  counter-clockwise  from  East 
(x  axis)  and  spans  0  to  360  degrees.  The  rings  in  the  plot 
represents  inclination  angles  (0)  ranging  from  0  (upward 
in  the  z  axis  direction)  to  90  degrees  (horizon). 

A  four-channel  RF/IF  front  end  is  used  to  collect  and 
digitize  the  outputs  from  each  antenna  element.  Figure  3 
shows  the  block  diagram  of  a  single  channel  of  the  RF/IF 
front  end.  The  digitized  outputs  from  each  channel  are 
stored  for  beam  forming  and  software  receiver  processing. 


GPS  signals  in  controlled  jamming  environment  were 
collected  in  an  anechoic  chamber.  Figure  1  shows  the 
setup  of  the  experiment.  A  two  by  two  conventional  GPS 
L1/L2  patch  antenna  array  was  placed  in  a  geodesic  dome 
as  shown.  The  antenna  elements  are  marked  with  number 
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adaptive  processor 


3.1.  MULTI-ANTENNA  ELEMENTS  GPS 
RECEIVER  INPUTS 


Fig.  2.  Sky  view  of  satellites  and  interference  sources  in 
the  experimental  setup 


Antenna 


Fig.  3.  Block  diagram  of  a  single  channel  GPS  receiver 
RF/IF  front  end 

3.  BEAM  FORMING  ALGORITHMS 

A  simple  spatial  adaptive  processor  as  shown  in 
Figure  4  is  used  in  this  study.  Inputs  from  the  antenna 
elements  k  are  applied  with  weight  (w^,  k 

=1,...K)  and  combined  to  generate  an  output  y. 

K 

(1) 

k=l 

where, 

w  =  W2  •••  =  ^2  •••  - 

Beam  forming  algorithms  are  designed  to  generate  a 
set  of  optimized  weight  w,  so  that  y  is  a  much  improved 
version  of  the  input  signal  x  in  that  the  interference  and 
jamming  signals  are  suppressed,  while  the  GPS  signals 
are  enhanced. 


The  signal  impinging  on  an  antenna  element  is  the 
collection  of  inputs  from  all  GPS  signals  in  direct  view  of 
the  antenna,  their  multipath  and  interference,  and  noise: 

P  Q 

Xk=Y^ gp^kp i^p 4p)+Y. (2) 

p=l  q=\ 

where, 

P\  Number  of  GPS  signals  received  by  the  antenna. 

Q\  Number  of  interference  (including  multipath)  signals 
received  by  the  antenna. 

Op,  (pp.  The  GPS  satellite  inclination  and  azimuth  angle. 

bkq;  The  k^^  antenna  spatial  phase  delay  factor  for  a 
signal  arrival  from  a  specified  direction. 

Vq'.  The  interference  signal  at  the  reference  location. 

Oq,  ^q\  The  interference  source  inclination  and  azimuth 
angle. 

Noise  component  for  the  antenna. 
gp{t)\  The  GPS  satellite  signal  received  at  a  reference 
location. 

00 

f^hpn)upt-nT^)  (3) 

n=-co 

In  Equation  (3),  fp-fn+f^p  is  the  GPS  signal 

carrier  frequency,  fn  is  GPS  satellite  LI  band  center 
frequency,  Ap,  (pp,  hp,  and  Up{t-nT^  are  the  satellite 
signal  Doppler  frequency,  amplitude,  carrier  phase, 
navigation  data,  and  CA  code  with  chip  duration 

Equation  (2)  can  be  written  in  a  more  compact  form: 
x  =  as +  bv  +  n  (4) 

where, 

■^  =  kl  gl  •••  gpJ’V=[vi  Vj  •••  Vq}  ’ 

n  =  [n^  •  •  •  ^kY 

a^\a^  •••  apYb  =  ^^  b^  ••• 

ap(p=l,  P)  and  b^(q=l,  ...,  Q)  are  the  spatial 

signatures  (also  called  steering  vector  or  directional 
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vector)  of  GPS  signal  source  p  and  interference  source  q. 
Assuming  all  signals  are  from  far  field  sources  and  that 
the  time  it  takes  for  a  signal  to  travel  between  array 
elements  is  much  smaller  that  inverse  of  the  receiver 
bandwidth,  the  spatial  signature  for  a  GPS  satellite  is: 


aip 

^j2^fpE\p 

II 

dip 

1 

^2- 

...  « 

_ 1 

^j^TrfpTKp 

(5) 


where,  t^=  g^-(r,-7^)l c,  g/is  the  unit  veetor 
pointing  from  the  signal  source  toward  the  A:*  antenna, 
and  are  position  vectors  of  the  A*  antenna  and  the 

reference  location,  respectively,  and  c  is  the  speed  of 
light.  For  the  4-element  antenna  shown  in  Figure  1,  if 
antenna  element  1  is  the  reference  location,  the  steering 
vector  associated  with  a  particular  signal  source p  is: 

"1 


a 


p 


jlTjfp  sin  6p  (cos^zi^  -sin^zi^  )d  ! c 


sin  6'  (cos^z^  +sin<z^  )f//c 


(6) 


e 


-j^¥p  sin 6'^  ^m(j)pdlc 


A  9-element  antenna  array  is  used  in  simulation 
studies.  The  9-element  antenna  array  assumed  the  layout 
as  shown  in  Figure  5.  The  steering  vector  for  this  antenna 
array  layout  is: 

sin  Op  (cos  (l)p  -sin  (l)p)dlc 


e 


-j'^Tfp  sin 6'^  ZQ^^pdlc 


e 


-jlTfp  ^m6p  (cos^zi^  +sin<zi^  )d ! c 


e 


jlTfp  sin 6'^  %m^pdlc 


1 


-j^Tf  sin 6'  sin^  Ic 


j27f  sin6'  (cos^zi  +sin^zi  )d ! c 


J'^^p  sin  6'  cos  (!)  die 


J'^4p  sin  6'  icoscj)  -smej)  )d  !  c 


(7) 


Fig.  5.  Layout  of  a  9-element  antenna  array  used  in 
simulation  study. 

3.2.  MOP  METHOD 

Assuming  that  GPS  signals,  interference,  and  noise 
are  uncorrelated,  the  total  combined  average  output  signal 
power  at  the  receiver  is: 

Ei\y{tf}  =  w"R^w  =  w"{R,  +  Ry  +  R„)w  (8) 

where,  R^=E{xx^},  R^=E{SS'^},  Ry  =E{VV"},  and 
=  E{nn^}  are  the  expected  correlation  matrix  for  the 

total  input  signal,  GPS  signals,  interference,  and  noise 
respectively.  In  the  MOP  method,  the  weight  vector  w  is 
obtained  by  minimizing  the  total  output  power  given  by 
(8),  while  constrain  the  gain  along  the  direction  of  known 
satellites: 

minw^R^w  subject  to  =  I  (9) 

The  solution  that  satisfies  (9)  is  (Godara,  2004): 

w  =  RMa^R;'ay  (10) 

where  a  contains  the  spatial  signature  for  those  satellites 
whose  directions  are  known  and  whose  signals  are  of 
interests  to  the  user.  Zoltowski  and  Gecan  (1995)  pointed 
out  that  the  constraints  applied  in  the  optimization  process 
consume  the  limited  number  of  degrees  of  freedom 
associated  with  a  phased  array.  It  may  provide  up  to  K-P 
spatial  nulls  to  cancel  interference.  As  can  be  seen  from 
the  analysis  results  to  be  presented  in  Section  5,  this  is 
truly  the  case  for  the  array  used  in  this  study. 

3.3.  SCORE  METHOD 

The  SCORE  method  was  proposed  by  Agee  et  al 
(1990)  and  was  applied  to  GPS  interference  cancellation 
by  Sun  and  Aimin  (to  appear).  The  SCORE  method 
explores  the  unique  characteristics  of  a  class  of  signals 
that  are  spectrally  self- coherent,  ie,  the  correlation 
between  a  signal  s{t)  and  its  frequency-shifted  version  for 
some  time  lag  is  non-zero: 

,  <s{t)s\t-z)e-^^’^^>^  7?f(r),Q(ll) 

f  KOr  >„<  \s\t-T)e-^^’'y  >„ 
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where  the  notation  (•)oo  represents  infinite  time  averaging 

operation.  The  function  Ps  (^)  is  referred  to  as  the 

spectral  self-coherence  function  of  s{t). 
7^^^(r)andi^^(0)  are  the  cyclic  autocorrelation  function 

and  autocorrelation  function  of  respectively.  The 
SCORE  algorithms  aims  at  obtaining  optimized  weight 
vector  w  that  maximize  the  cyclic  components  of  a 
receiver  input.  If  the  interference  and  noise  at  a  GPS 
receiver  are  not  spectrally  self-coherent  at  frequency 
separation  p,  the  cyclic  autocorrelation  of  the  receiver 
input  X  is: 

(r)  =  \a\Rf{T)  +  (r)  +  (r)  =  (r)  (12) 

The  CA  code  of  GPS  signals  Up{t-nT^  as  shown  in 
Equation  (3)  has  a  chip  duration  of  =  P77  ns  and  the 
code  period  T  is  7  ms.  The  GPS  navigation  data  has  a 
data  rate  of  5^  Hz.  As  a  result,  the  CA  code  repeats  itself 
20  times  within  one  navigation  data  bit.  Within  the  same 
navigation  data  bit,  a  GPS  satellite  signal  is  spectrally 
self-coherent  without  frequency  separation.  The 
collection  of  all  GPS  satellite  signals  at  a  receiver  input  is 
also  spectrally  self-coherent  without  frequency  separation. 
Because  of  this  property.  Sun  and  Aimin  (to  appear) 
developed  a  self-coherent  beam  former.  Two  set  of  inputs 
are  used  in  the  processor:  the  direct  input  x  from  the 
antenna  array,  and  a  delayed  version  of  the  input  x{t-mT) 
where  m  is  an  integer.  The  direct  input  x  is  applied  with  a 
set  of  weight  w  to  generate  an  output  y,  while  the  delayed 
input  x{t-mT)  is  applied  with  another  set  of  weight  w’  to 
generate  a  reference  output  y  \  The  difference  between  y 
and  y  e,  is  minimized  to  generate  optimized  weight  w 
This  minimization  is  done  in  the  least-square  sense.  The 
main  beam  former  generates  weight  w  by  maximizing  the 


correlation  between  y  and  y 

1 

^(14-,^')=  A’ 

(13) 

RyRy 

Ryy  =  E{yy')  =  w‘^R^{mT)w 

(14) 

Ry  =  E{yy}  =  w"i?^(0)w 

(15) 

R,  =  E{yy'}  =  w'^  RM^' 

(16) 

R^mT)  denotes  the  correlation  of  x  with  its  ml 
delayed  version  and  (0)  is  the  auto-correlation  function 

of  X.  Assuming  that  interference  and  noise  are 
uncorrelated  between  themselves  and  their  ml  delayed 
version,  then  RpmT)  =  R^-  Maximizing  F{w,  w ')  is 

therefore  the  maximization  of  total  GPS  signal  output. 

Substituting  (14),  (15),  and  (16)  into  (13)  and  solve 
for  the  optimization  problem,  the  desirable  weight  vector 
w  is  the  eigenvector  corresponding  to  the  largest 
eigenvalue  of  the  generalized  eigenvalue  problem: 

=  ^^^RAmT)R:^RymT)w  (17) 

MatLab  is  used  to  implement  both  the  MOP  and  the 
SCORE  algorithm.  The  correlation  matrix  R^  and  RxijnT) 


are  calculated  using  multiple  blocks  of  digitized  GPS 
input  data.  This  is  necessary  to  reduce  the  impact  of 
possible  navigation  data  bits  transition  within  an  input 
data  block  or  its  delay  version.  Figure  6  explains  the 
general  idea  of  data  block  selection.  Each  block  contains 
exactly  one  CA  code  period  of  data.  If  Data  Block  2 
which  contains  a  navigation  data  bit  transition  is  used  to 
calculate  the  correlation  functions,  the  self-coherence 
property  of  the  signal  will  be  non-existence  or  greatly 
reduced.  Using  the  averaged  correlation  results  of 
multiple  blocks  of  data  and  their  corresponding  delayed 
blocks  will  lessen  the  impact  of  the  navigation  data  bit 
transition. 


Navigation  data  bit  transition 


Data  Block  1 

Data  Block  2 

Data  Block  m 

Data  Block  m+1 
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Fig.  6.  Data  Block  Selection  for  SCORE  Algorithm 

4.  SIMULATION  RESULTS 

Simulation  inputs  containing  GPS  signals  and  noise 
only  are  generated  to  test  the  beam  forming  algorithm 
performance  in  this  project.  Inputs  from  two  different 
antenna  arrays  are  simulated.  Both  MOP  and  SCORE  are 
used  to  process  the  simulated  inputs  and  will  be  discussed 
in  the  following  two  subsections. 

4.1.  MOP  METHOD 

Figure  7  shows  two  example  antenna  gain  patterns 
generated  for  an  input  signal  that  contains  all  eight 
satellites  shown  in  Figure  2  using  the  MOP  method.  A  4- 
element  antenna  is  used  in  the  simulation.  A  single  beam 
constraint  is  applied  in  both  cases.  The  figure  shows  that 
MOP  is  capable  of  generating  a  well  defined  beam 
towards  the  given  satellite  direction. 

The  performance  of  the  MOP  method  degrades  as  the 
number  of  beam  constraints  increases.  This  is  evident  as 
shown  in  Figure  8  where  3 -beam  constraints  and  8-beam 
constraints  are  applied  for  the  4-element  antenna  input. 

The  number  of  antenna  elements  increases  the  degree 
of  freedom  of  the  MOP  beam  former.  Figure  9  is 
generated  using  the  same  input  as  that  of  Figure  8  and 
applied  to  a  9-element  planar-layout  antenna  array. 
Comparison  of  Figure  9  with  Figure  8  shows  clearly  the 
improvement  of  the  beam  forming  algorithm  performance 
due  to  the  increase  in  the  number  of  antennas. 
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Fig.  7.  Simulated  antenna  gain  pattern  using  4-element 
antenna  and  MOP  algorithm  for  single  satellite,  no 
jammer  scenarios. 


90 


Fig.  9.  Simulated  antenna  gain  pattern  using  9-element 
antenna  and  MOP  algorithm  for  multiple  satellites,  no 
jammer  scenarios. 

4.2.  SCORE  METHOD 

Figure  10  shows  the  antenna  gain  patterns  generated 
using  SCORE  algorithm.  The  input  signals  for  the  four 
patterns  contain  one,  two,  three  and  eight  GPS  signals 
respectively.  The  antenna  array  used  for  the  simulation 
has  4-elements.  The  SCORE  algorithm  generates 
satisfactory  results  for  one  and  two  satellite  scenarios. 
For  the  three  and  eight  satellite  scenarios,  the  algorithms 
can  not  operate  optimally. 

The  performance  of  the  SCORE  algorithm  does 
improve  as  the  number  of  antenna  elements  increases. 
Figure  11  is  the  result  of  using  a  9-element  antenna  array 
for  the  same  simulation  signals  used  to  generate  the  last 
two  patterns  shown  in  Figure  10. 

The  addition  of  jammers  will  further  influence  the 
beam  pattern.  Future  simulation  will  be  conducted  to 
study  this  problem. 


Fig.  8.  Simulated  antenna  gain  pattern  using  4-element 
antenna  and  MOP  algorithm  for  multiple  satellites,  no 
jammer  scenarios. 
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Fig.  10.  Simulated  antenna  gain  pattern  using  4-element 
antenna  and  SCORE  algorithm  for  single  and  multiple 
satellites,  no  jammer  scenarios. 


Fig.  11.  Simulated  antenna  gain  pattern  using  9-element 
antenna  and  SCORE  algorithm  for  multiple  satellites,  no 
jammer  scenarios. 

5.  EXPERIMENTAL  DATA  ANALYSIS 

A  total  of  33  sets  of  experimental  data  containing 
GPS  and  jamming  signals  were  taken  for  this  study.  The 
following  subsections  will  present  the  results  of  the  MOP 
and  SCORE  beam  forming  algorithms  and  software 
receiver  acquisition  performance  based  on  the 
experimental  data  processing. 

5.1.  MOP  METHOD 

Among  the  33  sets  of  experimental  data,  there  are  6 
sets  containing  single  FM  chirp,  6  sets  containing  single 
BPSK,  and  21  sets  containing  single  and  multiple 
broadband  interference  sources.  Figures  12  through  20 
summarize  the  beam  forming  and  receiver  acquisition 
results  of  these  experimental  data.  Figure  12,  14,  16,  and 
1 8  compares  the  average  number  of  successful  GPS  signal 
acquisition  in  the  presence  of  FM  chirp,  BPSK,  single 
broadband,  and  multiple  broadband  interference  sources, 
while  Figure  13,  15,  17,  and  19  compares  the  post¬ 
acquisition  signal  to  noise  ratio  for  the  successfully 
acquired  signals.  It  is  evident  from  these  figures  that  the 
MOP  method  worked  effectively  to  cancel  FM  chirp  and 
BPSK  interferences  at  all  interference  power  level  tested. 
The  method  is  also  functional  to  some  degrees  when 
dealing  with  single  broadband  jamming  source.  For 
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multiple  broadband  jamming  sources,  its  performance 
degrades  considerable  at  the  30  and  40  dB  JSR  levels. 


Comparison  of  Average  Successful  Acquisition  for  FM  chirp 
Interference  (Total  8  Satellites  Available) 


Fig.  12.  Average  number  of  successful  GPS  signal 
acquisition  in  the  presence  of  FM  chirp  interference 


Comparison  of  Post  Acquisition  Average  SNR  for  FM  Chirp 
Interference 
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Fig.  13.  Average  post-acquisition  GPS  signal  to  noise 
ratio  in  the  presence  of  FM  chirp  interference 


Comparison  of  Average  Successful  Acquisition  for  BPSK 
Interferences  (Total  8  Satellites  Available) 
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Fig.  14.  Average  number  of  successful  GPS  signal 
acquisition  in  the  presence  of  BPSK  interference. 


Comparison  of  Post  Acquisition  Average  SNR  for  BPSK  Interference 
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Fig.  15.  Average  post-acquisition  GPS  signal  to  noise 
ratio  in  the  presence  of  BPSK  interference 


Comparison  of  Successful  Acquisition  for  Broadband  Interference 
(Average  Single  Jammer) 
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□  Ave.  Single  Channel  Output 
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Fig.  16.  Average  number  of  successful  GPS  signal 
acquisition  in  the  presence  of  single  broadband 
interference. 


Comparison  of  Post  Acquisition  Average  SNR  for 
Broadband  interference  (Singie  Jammer) 
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Fig.  17.  Average  post-acquisition  GPS  signal  to  noise 
ratio  in  the  presence  of  single  broadband  interference 
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Comparison  of  Successful  Acquisition  for  Broadband  Interference 
(Average  Multiple  Jammer) 


that  using  single  beam  constraint,  while  the  three-beam 
constraints  cannot  generates  satisfactory  weighted  output. 
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□  Ave.  Single  Channel  Output 


■  Beam  Forming  Output 
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Fig.  18.  Average  number  of  successful  GPS  signal 
acquisition  in  the  presence  of  multiple  broadband 
interferences. 
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Comparison  of  Post  Acquisition  Average  SNR  for 
Broadband  interference  (Muitipie  Jammers) 
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Fig.  19.  Average  post-acquisition  GPS  signal  to  noise 
ratio  in  the  presence  of  multiple  broadband  interferences. 


File  160615  SUV  10  SNR=25.32 


(a) 


The  following  groups  of  figures  (Figure  20  (a)-(f)) 
show  some  examples  of  the  antenna  gain  pattern  and 
receiver  acquisition  results.  Antenna  gain  pattern  and 
software  GPS  receiver  acquisition  stage  correlation  code 
phase  output  for  experimental  data  collected  using  a  4- 
element  antenna  are  plotted  in  this  group  of  figures. 
Figures  20  (a),  (b),  and  (c)  show  the  results  for  single 
jammer,  single  beam  constraint  for  40  JSB  FM  chirp, 
BPSK,  and  broadband  respectively.  Figure  20  (d)  and  (e) 
shows  the  outputs  involving  two  and  three  broadband 
jammers.  Broadband  jammer  sources  with  JSR=40  dB 
are  included  in  the  data  for  both  figures.  A  single  satellite 
is  used  to  constraint  the  beam  direction.  In  Figure  20  (f), 
a  single  broadband  jammer  (JST=40  dB)  is  included  in 
the  data  and  two  satellites  are  used  to  constraint  the  beam 
forming  algorithms. 

Additional  processing  is  performed  for  two  and  three 
beam  constraints  using  the  MOP  method.  The  general 
conclusion  is  that  post-acquisition  signal  to  noise  ratio  of 
satellites  using  two-beam  constraints  is  compatible  with 
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Fig.  20  (a)-(f).  Antenna  gain  patterns  and  receiver 
acquisition  results  for  experimental  data  collected  using  a 
4-element  antenna  array.  The  results  is  obtained  using 
MOP  method. 


5.2.  SCORE  METHOD 

The  SCORE  method  did  not  generate  outputs  that 
lead  to  consistent  acquisition  results.  Figure  21  shows 
four  sample  antenna  gain  pattern  generated  for  data 
containing  a  variety  of  jamming  sources.  Figure  21(a), 
(b),  (c),  and  (d)  each  contains  a  single  jammer  with 
JSR=20  dB,  a  single  jammer  with  JSR=40  dB,  two 
jammers  with  JSR=40  dB  each,  and  three  jammers  with 


JSR=30  dB  each,  respectively.  Broadband  jammers  are 
contained  in  all  four  data  files.  As  can  be  seen  that  the 
SCORE  beam  forming  algorithm  places  nulls  in  an 
inconsistent  manner  in  the  patterns.  This  is 
understandable  since  there  are  a  total  of  eight  satellites  in 
the  experimental  setups  while  only  4  antenna  elements  are 
used  to  collect  the  data. 
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Fig.  21  (a)-(d).  Antenna  gain  pattern  for  experimental 
data  collected  using  a  4-element  antenna.  SCORE 
algorithm  is  used  to  obtain  the  results. 

6.  CONCLUSIONS  AND  FUTURE  WORKS 

This  paper  presented  results  of  integrating  two  digital 
beam  forming  algorithms,  minimum  output  power  and 
SCORE,  with  a  software  GPS  receiver  to  study  the 
effectiveness  of  the  beam  forming  algorithms.  Both 
experimental  and  simulation  data  are  used  to  test  and 
evaluate  the  performance  of  these  algorithms.  The 
software  GPS  receiver  is  used  to  perform  signal 
acquisition  of  the  beam  former  output.  Our  study  showed 
that  MOP  can  be  an  effective  interference  cancellation 
technique  when  satellite  signal  angle  of  arrival  is  known 
and  when  enough  degrees  of  freedom  are  available.  Our 
experimental  results  show  that  the  MOP  method  works 
with  all  3  types  of  interference  sources  with  JSB  up  to  40 
dB. 

Our  experience  with  the  SCORE  method  suggests 
that  it  requires  large  numbers  of  antenna  elements  when 
the  number  of  available  satellites  is  large.  It  will  produce 
erroneous  results  if  a  jammer  also  has  spectrally  self- 
coherent  properties.  Because  of  its  sensitivity  to  the 
special  spectral  self-coherent  properties  of  the  CA  code,  it 
may  be  used  for  ground-based  jammer  identifications. 

The  study  shows  that  integrating  beam  forming 
algorithm  with  Software  GPS  receiver  can  provide  a 
powerful  means  to  develop  anti-jam  receivers. 

For  future  works,  more  experiments  involving 
different  types  of  antenna  elements,  antenna  layout 
configurations,  more  receiver  input  channels,  jammer  and 
satellite  configurations,  jammer  types  and  power  levels 
should  be  carried  out.  Additional  simulation  works  need 
to  including  jammers.  Analysis  using  tracking  and 
position  calculations  will  also  be  investigated  in  the 
future. 
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